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1. Introduction 

Recently, a lot of new experimental reports are delivered on charmed baryons [|]]. BaBar 
and Belle give very accurate results, and precision analysis can be accomplished. In addition, new 
experiments, such as J-PARC, PANDA, LHCb, and Belle II are coming and expected to give further 
insight into charmed baryons. 

Mass spectrum of the singly charmed baryons is determined in the high precision by experi- 
ments. Experimental status of the ground state is three or four-star, evaluated by the particle data 
group. The excited states are also investigated fairly well. 

In contrast to the singly charmed baiyons, experimental data for the doubly and triply charmed 
baiyons are not sufficient. Only one candidate for the doubly charmed baryon. Sec, has been 
reported by the SELEX Collaboration [^, while Sec is not confirmed in the other experiments 
such as BaBar ^ and Belle [Q] groups. Experimental and theoretical cross checks are needed to 
establish Sec. The other doubly charmed baryons and the triply charmed baryon have not been 
found by experiment, yet. Theoretical predictions for the doubly and triply charmed baryons are 
helpful to discover these states. 

So far, almost all lattice full QCD calculations for charmed baryon spectrum have been per- 
formed with the staggered dynamical quarks ^, ^ ^. They use some technique, such as con- 
verting the staggered propagators or mixed action, for the valence light quarks to deal with the 
tangled flavor structure of the staggered quarks. Although these approaches have a correct contin- 
uum limit, it breaks the unitarity and complicates the continuum extrapolation. It is preferred to 
use the other type of dynamical quarks which is simple in flavor. Another point is that their chiral 
extrapolations suffer from large higher order corrections. The pion masses are limited to 220 - 290 
MeV. NLO SU(2) heavy baryon chiral perturbation theory is employed to extrapolate their data 
to the physical point, but it shows a bad convergence even with m^j = 220 MeV. It is desirable to 
perform a simulation directly on the physical point. 

The ETMC group studied the charmed baryons with Nf = 2 dynamical twisted mass quarks 
and Osterwalder-Seiler strange and charm valence quarks [|l^]. They found S^c = 3.513(23)(14) 
GeV, which agrees with the SELEX experimental value E^^r^^^ = 3.519(1) GeV. This is the only 
result that is consistent with the SELEX experiment. The other lattice QCD simulations show 
deviations from it. This disagreement in lattice QCD must be resolved. A subtle issue in the 
ETMC calculation is that the heavy quark mass conection may not be under control at their lattice 
spacings, a = 0.09 — 0.06 fm. Their results for charmed baryons, especially for the triply charmed 
baryon Q.ccc, do not show a cleai^ scaling behavior. It is better to employ a heavy quai^k action that 
handles mass dependent lattice artifacts in the formulation. Chiral extrapolation of ETMC data 
from nijc = 260 MeV using the NLO heavy baiyon chiral perturbation theory is also problematic. 

In Ref. [11], we have shown that the charm quark mass corrections are under control at = 
2.194(10) GeV by adopting the relativistic heavy quark action of Ref. [[T^. It removes the leading 
cutoff errors of 0{{mQa)") and the next to leading effects of 0{{mQa)" {aAgcD)) for arbitrary order 
n. We calculated the spectra of mesons involving charm quarks using the relativistic heavy quark 
action on the 2-1-1 dynamical flavor PACS-CS configurations of 32-^ x 64 lattice reweighted to the 
physical point for up, down and strange quark masses. We found our results are consistent with 
experiment at a percent level, and so are those for the decay constants with a few percent accuracy. 
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1.90 
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0.13663375 


80 


2000 



Table 1: Simulation parameters. MD time is the number of trajectories multiplied by the trajectory length. 



Based on this result, we extend our calculations to the charmed baryon sector. The notable 
point is that our measurements are performed at the physical point. We are free from the conver- 
gence problem of the heavy baryon chiral perturbation theory. We check if our masses of singly 
charmed baryons reproduce the experimental values to confirm validity of our calculation. We also 
evaluate the doubly and triply charmed baryon spectrum as our predictions. 



2. Setup 

Our calculation is based on a set of 2 + 1 flavor dynamical lattice QCD configurations gener- 
ated by the PACS-CS Collaboration on a 32^ x 64 lattice using the nonperturbatively 0{a)- 
improved Wilson quark action with = 1.715 [14] and the Iwasaki gauge action at /3 = 1.90. 



The aggregate of 2000 MD time units were generated at the hopping parameter given by (k"''^, k'^) = 
(0.13778500,0.13660000), and 80 configurations at every 25 MD time units were used for mea- 
surements. We then reweight those configurations to the physical point given by (K"„t/,K;) = 
(0.13779625,0.13663375). The reweighting shifts the masses of 7i and K mesons from rtij^ = 
152(6) MeV and niK = 509(2) MeV to = 135(6) MeV and uik = 498(2) MeV, with the cutoff 
at the physical point estimated to he. a^^ = 2.194(10) GeV. Our parameters and statistics at the 
physical point are collected in Table |I[ 

The relativistic heavy quai^k formalism [ p^ is designed to reduce cutoff errors of 0((mQa)") 
with arbitrary order n to 0{f{mQa){aAQCD)^), once all of the parameters in the relativistic heavy 
quark action are determined nonperturbatively, where /{mga) is an analytic function around the 
massless point mga = 0. The action is given by 



Sq = Y,QxD.,yQy, 

x,y 

D,,y = 5^. - KqY, \{rs - vri)U,j5^^.y + (r, + VT^-)^),,-5,,y+? 



(2.1) 



(1 - 74)t/.v,45,+4,^, + (1 + 74)t/j45,,,+4 



22 Fij {x) Oij +CeY, ^'-4 {x) <^i4 



'J 



(2.2) 



The parameters rs,CB,CE and v have been adjusted in Ref. Qllp. It should be noticed that the 
pai^ameter v is determined non-perturbatively to reproduce the relativistic dispersion relation for 
the spin-averaged 15 state of the charmonium. The heavy quark hopping parameter Kq is set to 
reproduce the experimental value of the mass for the spin-averaged IS state. Our parameters for 
the relativistic heavy quark action are summarized in Table |^ 

We use the relativistic operators to obtain charmed baryon spectrum, because the relativistic 
heavy quark action is employed in our calculation. Charmed baryons can be classified under 4 x 4 x 
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rs Cb 


Ce 


0.10959947 


1.1450511 


1.1881607 1.9849139 


1.7819512 



Table 2: Parameters for the relativistic heavy quark action. 



4 = 20 + 20i + 2O2 + 4. In addition to / = 3/2 decuplet-type 20-plet, there are 7 = 1 /2 octet-type 
20-plet and 4-plet. 

7=1/2 octet-type baryon operators are given by 

0ff{x) = e"'^i{q}{x)fCYsq',{x))qUx), (2.3) 
C = 7472, (2.4) 

where f,g,h are quark flavors and a,b,c are quark colors, a = 1,2 labels the z-component of the 
spin. The £-type and A-type are distinguished as 

Z-type : — , (2.5) 

A-type : — , (2.6) 

where o'-^^l'' = O^^^'' - 0^f'\ 

The decuplet-type J = 3/2 baryon operators are expressed as, 

D^fjx) = e"'%{q}{x)fCr^q'i{xM,{x), (2.7) 



= e"'^[{{q}{ 



D{fjx) = e^''^[i{q}ix)fcroq'jxMAx) 



-{{q}{x)fcr^q'^{x))ql,{x)]/3, (2.8) 
D{f/,(x) = e"'^[{{q}ix)fCroq'^{x))q^A^) 

-{{q}{x)Ycr^q'^{xMAx)]/3, (2.9) 

Z){f)2W = e"''i{q%x)fcr^q'^^ix))q',,ix), (2.10) 



r± = (7iT/72)/2,ro = 73. (2.11) 



The baryon correlators are calculated with exponentially smeared sources and a local sink. 
The smearing function is given by ^{r) = Aexp(— Br) at r 7^ and ^{0) = 1. We set A = 1.2, 
B = 0.07 for the ud quark, A = 1.2, B = 0.18 for the strange quark, and A = 1.2, B = 0.55 for 
the charm quark. The number of source points is octupled and polarization states are averaged to 
reduce statistical fluctuations. Statistical eiTors are analyzed by the jackknife method with a bin 
size of 100 MD time units (4 configurations), as in the light quark sector [|l^]. We extract charmed 
baryon masses by fitting correlators with exponential functions. 

3. Singly charmed baryon spectrum 

Our results for the singly charmed baryon spectrum at the physical point are summarized 
in Fig. [1} All our values for the charmed baryon masses are predictions, because the physical 
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Figure 1: Our results for the singly charmed baryon spectrum (left panel), and comparison of mass with 
other lattice QCD results (right panel). 



chaim quark mass has already been fixed with the charmonium spectrum. We found the predicted 
spectrum is in reasonable agreement with experiment. We also compare our value for A^. with other 
results by recent lattice QCD simulations using the dynamical staggered quarks [Q, ^ ^, and the 
twisted mass quarks |[T0|]. All results are consistent with each other. 

Fig. ^ displays several mass differences. We have consistent results with experiments in 2 a 
accuracy. The decomposition of 7 = 1/2 Z-type and A-type baryons, as well as that of 7 = 3/2 and 
7=1/2 charmed baryons, are successful. 

It is noted that several systematic errors have not been evaluated, yet. One is finite size effects. 
Though NLO heavy baryon chiral perturbation theory predicts that finite size effects for charmed 
baryons are less than 1 %, higher order terms can give significant contributions. A direct lattice 
QCD check by comparing spectrum on different lattice volumes is desirable. Another aspect is that 
strong decays such as Ec — )■ Ac?! are not taken into account. Sc — ^ AcTT is kinematically prohibited 
on our lattice. Our estimates should be considered as the lower mass limits for unstable baryons. 
Moreover, we have not performed the continuum extrapolation. A naive order counting implies 
that the cutoff effects of O ( f{mQa) {oAqcd ) , /("iga) {a^QCD ) ^ ) from the relativistic heavy quark 
action are at a percent level. Additional calculations are needed to remove these systematic errors. 

4. Doubly and triply charmed baryon spectrum 

For doubly and triply charmed baryons, an experimental value has been reported only for 
Sec, although the experimental status is controversial. In the other channels, lattice QCD gives 
predictions to experiments. 

Fig. ^ shows our results for the doubly charmed baryons. Our estimate for ms„, clearly deviates 
from the experimental value by SELEX Collaboration The difference is 4a. We compare our 
result for m^^,^ with those by other lattice QCD calculations. We have a consistent value with other 
lattice QCD calculations, except for that by ETMC. 

Similarly, Fig. ^ displays lattice QCD results for the triply charmed baryon from several 
groups. Our prediction agrees with those by others, except for ETMC value. We also plot niQ,^^^ 
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Figure 2: Comparison of mass differences of 
panel), E* — types (lower panel). 



Ac types (upper left panel), E* — types (upper right 



— 3/2 mj/^ mass difference. A slight discrepancy is observed. For a more definite comparison, 
evaluation of systematic errors is necessary. We need to take the continuum extrapolation. 
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